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ABSTRACT 

The spacing of the acoustic peaks in the cosmic microwave background radiation 
anisotropy multipole spectrum has been claimed to provide the value of the total cos- 

Q^ I mological density overtly, "written on the sky." Through a semianalytic analysis of the 

cosmological evolution of the sound horizon and the physics of decoupling we address the 
robustness of the relation between the peak spacing and the cosmological density. In fact, 

Q . the asymptotic distance and horizon scalings often used are not good approximations, 

and the individual densities and equations of state of different components do enter the 
problem. An observed spacing could be fit by models with different total densities. We 
investigate the different regions of density-equation of state parameter space and also pro- 
^ ■ vide accurate fitting formulas for the peak spacing as a function of matter density, total 

0> ! density, and additional component equation of state (e.g. cosmological constant or cosmic 

strings) . Limits provided by peak spacing measurements on the number of neutrino species 

CN ' and the baryon-photon ratio are also addressed. 

^H ' 

0> ■ 1. Introduction 

p . The cosmic microwave background (CMB) radiation provides a cornucopia of cosmo- 
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logical information about the early universe and the formation of large scale structure. 



H i With the increasing sensitivity of medium scale angular anisotropy experiments, both 

^ ■ ground based and the near future space missions MAP and Planck, further information 

on the properties of our universe is imminent. Measurements on scales of 10-60 arcmin- 
utes cover the horizon scales around the epochs of recombination and last scattering and 
r\ \ promise a wealth of data on the large scale properties of our universe as well as indications 

c^ I of the early universe origin of fluctuations that become cosmological structure. 

In particular, they hold the hope of measuring the total cosmological energy density 
and hence elucidating the ultimate fate of our universe, as well as providing evidence on the 
validity of the inflationary picture of perturbation generation. Such a lofty goal would be 
reached through a series of measurements using straightforward, well understood, physical 
theory, and would avoid more tortuous and uncertain paths such as the distance ladder 
and the nonlinear evolution of complex galaxies and stars. 

With the prospect of the revelatory data looming, it behooves us to analyze carefully 
the extent to which the translation from peak multipole spacing measurements to the sin- 
gle number of the total energy density of the universe is indeed obvious and clean. At 
first sight it seems astonishing that despite having different species of cosmological inhab- 
itants - photons, baryons, neutrinos, cold dark matter, cosmological constant, etc. - and 
an uncertain expansion rate (Hubble constant) that the spacing should depend on this one 
single parameter. Fortunately, the physics of baryon-photon decoupling and cosmological 
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distances is sufficiently simple that we can go a long way through analytic investigation. 
Of course, numerical calculations of not only the peak spacing but their amplitudes and 
the full multipole spectrum of microwave background anisotropies have been carried out 
(e.g. Zaldarriaga, Seljak, & Bertschinger 1997; Kamionkowski, Spergel, & Sugiyama 1994), 
and a rigorous derivation of cosmological parameter values should rely on the full informa- 
tion and accuracy they possess (e.g. Zaldarriaga, Spergel, & Seljak 1997; Bond, Efstathiou, 
& Tegmark 1997). But for understanding the physical origin of the results in different re- 
gions of parameter space, this analytic exploration provides a heuristic and surprisingly 
accurate portrait. 

2. Acoustic Oscillation Anisotropies 

On medium angular scales, 9 ~ 10' — 2° or multipoles / ~ 100 — 1000, the dominant 
anisotropy in the CMB is expected to be acoustic, or Doppler, oscillations (Hu, Sugiyama, 
& Silk 1997; Bond et al. 1994; occasionally these are called Sakharov (1966) oscillations). 
These are due to harmonic density perturbations in the coupled baryon-photon fluid in the 
prerecombination epoch. While after recombination these will grow to form the large scale 
structure of the universe, the density perturbations cannot grow as long as the matter is 
tightly coupled to the radiation, e.g. through Thomson scattering. This situation persists 
even after the end of the radiation dominated era due to the high entropy, or low baryon 
to photon ratio. Thus even when the matter energy density governs the expansion rate 
of the universe, the coupling forces the behavior to be oscillational - a sound or acoustic 
wave. 

As the universe expands, the Thomson scattering rate drops below the cosmological 
expansion rate and the coupling reactions become inefficient. Matter recombines to neutral 
atoms and the ionization fraction freezes out. The photons become free streaming and the 
lack of further interaction preserves the density irregularities imprinted on the photon 
field by the matter oscillations. These appear as angular anisotropies, with the largest 
angular scale arising from the largest matter wavelength A - that of the sound horizon at 
decoupling. 

In CMB studies one often works in multipole space / ~ 9~^ ~ A~^. A length scale 
translates to an angle 6*, or multipole /, through the angular diameter distance Va by 
/ = kva = 2nra/X = tt/^, where k is the wavenumber. Thus one expects a harmonic series 
of anisotropies corresponding to acoustic normal modes, with peaks aX li : I2 : h ■ ■ ■ ^ I : 
2:3.... The peak spacing A/ is a direct measure of the maximum wavelength, the horizon 
diameter A = 2r^. Thus, 

A/ = Ttra/rh. (1) 

In attempting to learn about cosmological parameters by observing the anisotropy 
acoustic peak spacing or locations, we are simply applying the classical angular diameter 
distance test to the CMB. This has many virtues over its use with galaxies and other 
conventional objects, including astrophysical simplicity and physical understanding of the 
source, lack of evolutionary or selection effects, and the high redshift improving the dis- 
crimination between cosmological models, as advocated by Linder (1988a). 

In order to compare observations of the multipole spacing to cosmological models 
we need three quantities: the angular diameter distance ra out to redshift z, the sound 



horizon scale r^ at redshift 2;, and the redshift of decouphng Zdec^ at which we evaluate 
these functions. Each depends on numerous subsidiary cosmological variables. Before 
looking at these dependences let us note a few issues regarding the applicability of this 
acoustic oscillation picture. 

The ansatz relies on the presence of harmonic perturbations in the coupled baryon- 
photon fluid prior to decoupling. While this is a generic feature of the adiabatic fluctuations 
generated in inflationary scenarios, it does not have to exist generally. The presence of 
isocurvature perturbations, such as those arising from topological defect models, suppresses 
some of the peaks and changes the peak spacing (see Hu & White 1996) . It is also possible 
to arrange active, causal processes in such a way as to create peaks that mimic the results 
of inflation (Turok 1996). 

Moreover, although in a purely harmonic mode the first peak occurs at multipole 
h = A/ and the spacing Im+i — Im = ^l for all m, not just to :» 1, gravitational forcing 
terms in the oscillations slightly change these predictions (Hu & White 1996). Given a 
theory for the origin of the perturbations, however, the locations and small to spacings 
(they reach the asymptotic value by m ~ 3) can be related to the pure harmonic value A/. 
Thus, though we phrase all our results in terms of the spacing A/, they can be applied as 
well to the peak locations. 

Although adiabatic perturbations arise naturally within inflationary theory, we will 
not restrict ourselves to flat cosmological models, both because of the possibility of open 
inflationary universes, and because adiabatic perturbations can be generated via other, 
causal processes. With these points in mind, this paper operates within the picture of the 
acoustic peaks directly tracing the sound horizon scale, for however far that is valid. 

3. Physical Variables 

As illustrated in (1), the peak multipole spacing depends on two distance scales and the 
redshift at which they are evaluated. The distance-redshift relations are straightforward 
cosmological expressions. The redshift is that of decoupling, for this is when the physical 
imprinting of the acoustic anisotropies in the CMB temperature pattern occurs, when the 
photons become unaffected by further interactions with the matter. 

3a. Angular Diameter Distance 

In a homogeneous and isotropic universe, i.e. a Friedmann-Robertson- Walker (FRW) 
model, the angular diameter distance is given by 



r^[z) = (l + z)-^(l-OT)"^/^sinh 



{l-nrY'^ p^ dy/H{y) 



(2) 



where Q.t is the total cosmological energy density in units of the critical density pc 
3iy^(0)/87r, y = 1 + z, and H{y) is the Hubble parameter. 



H{y) = HiO) 



J2^.y'^'+''^ + {I - nT)y' 
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(3) 



(See Linder 1988a,b, 1997 for general distance-redshift and cosmological parameter ex- 
pressions for multiple components and equations of state.) The equation of state of a 
cosmological component is given by cr = p/ p, where p is its pressure, p its energy density, 
and Oct its dimensionless density p/ Pc- 

Note that Va depends not only on Ot = XI ^cr but also on each component individually, 
though generally one particular equation of state will dominate at a given redshift. One 
subtlety involves the use of the FRW model despite the dumpiness of matter - i.e. structure 
- along the line of sight, at least at low redshift. The question of when to use clumpy model 
distances instead of FRW distances is not wholly settled, but according to recent criteria 
of Linder (1998), FRW relations can be applied with confidence when considering angular 
scales 9 > 5" (/ < 10^), as we do. 

3b. Sound Horizon Scale 

The particle horizon distance is r^ = a / dt/a, where a{t) is the expansion parameter, 
or scale factor, and t is the time. To obtain the sound horizon we simply correct this for 
the speed of propagation, or sound speed Cg, of the acoustic waves: 

/"OO 

th{z) = {1 + z)-^ dycM/H{y), (4) 

Jl + Z 

since H{y) = a~^da/dt and y = a~^. 

The sound speed is closely related to the equation of state parameter a: Cs = \fdipjdp 
while a = p/p. At the epoch of decoupling and earlier we expect the only two significant 
components to be nonrelativistic (pressureless) matter, a ~ {l/3)v'^/c^ ~ 0, and relativistic 
particles - photons and neutrinos - with a = 1/3. The only acoustic coupling that exists, 
however, is between the baryons and photons. Thus the total pressure fluctuation in 
the acoustic medium is dp = dpj = {l/3)dpj and the total energy density fluctuation is 
dp = dp^ + dpb- 

For adiabatic perturbations the total entropy fluctuation must vanish; since the en- 
tropy is directly proportional to the baryon-photon number ratio rj = Uh/n^, this im- 
plies {dn/n)b = {dn/n)-y. For matter, {dn/n)b = {dp/p)b while for blackbody radiation 
the number density goes as the temperature T^ and the energy density goes as T^, so 
(Sn/n)^ — (3/4) (5p/p)^. Thus, dpb = {3/4){pb/P'y) dp^. Finally, the evolution with expan- 
sion of the two components is pb ~ y^ and pj ^^ y^, so we flnd 



Deflning a redshift of baryon-photon equality, yeqib,-/) = ^b/^-y, we see that for y ^ 
^69(6,7)5 when photons dominate, the effective equation of state is aeff = c^ = 1/3, while as 
the universe expands and baryons come to dominate for y <^ yeqib,^)^ (^eff -^ 0. Thus the 
sound speed varies betweens l/v3 and 0, decreasing as baryons become more important 
(this can be viewed as increasing the effective mass of the acoustic phonon - see Hu, 
Sugiyama, & Silk 1997). Note that the redshift of baryon-photon equality is not the same 



as the redshift of matter-radiation equality, unless all matter is baryonic (no cold dark 
matter, for example), and all radiation is photons (e.g. no light or massless neutrinos). 

3c. Decoupling Redshift 

The last ingredient necessary is the redshift at which to evaluate the two distance 
scales, that of decoupling when the photons receive their final anisotropy imprint. The 
coupling process is dominated by Thomson scattering between the photons and the free 
electrons of the perennially ionized hydrogen and proceeds at a time dependent reaction 
rate T{z) = arneiz), where ar is the Thomson cross section and Ue the electron number 
density. 

As the universe expands, the photon gas cools, eventually dropping below the ion- 
ization energy of hydrogen, leading to recombination of the electrons into atoms. Due to 
the high entropy, or low baryon-photon number ry, this actually occurs when the photon 
temperature is 2 Inry ^ 1/40 of the hydrogen ground state ionization energy. However, al- 
though related, the decoupling of the photons and baryons is a slightly different situation: 
what is important here is the number of interactions ongoing between the two components. 
When the Thomson reaction rate drops below the expansion rate H ~ t~^ of the universe, 
or equivalently the photon mean free path exceeds the light horizon, then interactions 
effectively cease: Nint = L Tdt <1. 

The criteria for decoupling, therefore, is T < H (cf. Kolb & Turner 1990). Both rates 
are redshift dependent: as the universe expands (redshift decreases) the electron density 
is diluted and also the Hubble parameter decreases. Setting T = H defines the desired 
redshift of decoupling, Zdec- To find the electron density rig = Xe'r]n^{T^) in F one cannot 
simply use the Saha expression for the ionization fraction X^ = ne/ni, because the atomic 
levels are not in statistical equilibrium. Jones & Wyse (1985) provide an excellent analysis 
of the relevant physics and applicable approximations. We follow their work, generalizing 
it slightly. By comparing various atomic excitation rates, they reduce the expression 
for the ionization fraction as a function of redshift to a Riccati equation and find that 
near decoupling a quasiequilibrium WKBJ solution obtains, similar to the analysis for the 
nucleosynthesis rate equations done by Esmailzadeh, Starkman, & Dimopoulos (1991). 

This fixed point solution gives X^ ~ (Ot^^)^^^/~^(-s, ^t) (^b^^)~^, where h = 
H{0)/{100 km s~^ Mpc~^) and f{z, O^) involves the ratio of the decoupling redshift to the 
redshift of equality of energy density in the matter and radiation, Zeq- Retracing their anal- 
ysis for a general, not just dust (zero pressure), universe, i.e. allowing arbitrary equation 
of state components with density Qcr, reveals that the factor (OT/i^)^^^/~^(-s, f^r) arises 
from {dt/dz)~^ ~ H{z) and the Qbh^ factor comes from rib, since the total recombination 
rate is proportional to UeUp ~ nf while the competing total ionization rate is proportional 
to Ub (ne. Up are the number densities of electrons and protons, respectively). Therefore 
the generalized decoupling condition F = if is 

^{Zdec) = O-TUe ~ X^Ub ~ H{z)n^^nb = H{Zdec) ,^.. 

(6) 



Zdec Z/l^cyl^ J 



i.e. Zdec is practically a constant (except that the atomic transitions also depend on the 
background temperature T^). 



One has the remarkable conclusion that cancellations in H{z) and ni, lead to the 
redshift of decoupling being extraordinarily insensitive to the baryon content of the uni- 
verse, the total density and density of individual components, and the overall expansion 
rate. The redshift is almost completely determined by the photon temperature, which is of 
course fixed by CMB measurements, and implies Zdec — 1100. This is the value at which 
we evaluate the distances Ta and rh for all models. 

Jones & Wyse (1985) carefully examine the approximations they use and these appear 
robust. We will note that the leading correction term (the last term in their equation A13) 
is proportional to Q^/^t and so the approximation should be especially good for universes 
with Qh <^ Qt (as they state; also see underived equation 3.101 of Kolb & Turner 1990) - 
this holds for all models we consider. 

3d. Summary of Cosmological Dependences 

Before calculating ra{zdec) and rh{zdec) to obtain the multipole peak spacing we pause 
to review the cosmological ingredients that enter into the problem. In broad terms: r^ 
involves Ot and H{z); r^ involves Cs{z) and H{z); Zdec involves T{z) and H{z). Each of 
these major variables, however, has subsidiary elements entering: 

• H{z) - component densities Vta and equations of state a; current value h = 
if (0)/(100 km s-i Mpc-i) 

• Cs{z) - baryon-photon ratio 77; helium abundance Y through translation of ni, to 
Ofe (i.e. average mass of a baryonic particle: see Section 5b) 

• O1/3 - number of neutrino species N^,; Hubble constant h affects translation of T^ 
to O1/3 ~ h~'^T^] ^1/3 enters into H{z) and definition of Zeq = O1/3/O0 

• ra{z) - dumpiness of density distribution, as discussed in Section 3a 

We keep these in mind for two reasons: 1) to understand the original claims that 
the peak spacing A/ was degenerate in these ingredients and depended only on Oy, and 
2) to attempt to use the wealth and accuracy of CMB anisotropy data to probe these 
additional cosmological quantities when we see how the putative degeneracy breaks down. 
After obtaining results for A/ in the next section we return to this question in Section 5. 

4. Peak Spacing 

Using equations (l)-(5) and Zdec = 1100 we obtain analytic or quadrature expressions 
for the anisotropy multipole peak spacing A/. We examine open and flat models, Oj- < 1, 
with three cases of components: pure dust (cr = 0) with density Oq; dust plus a cosmological 
constant (cr = —1) with densities Oq, 0_i; dust plus a cr = —1/3 component with densities 
Oo, 0_i/3. [0_i/3 does not include curvature energy 1 — Oy; it is a distinct component 
with p = — (l/3)p, corresponding perhaps to a cosmic string network. Its function is simply 
to illustrate the effect of the equation of state on the anisotropy multipole pattern. It is 
also interesting to consider because it doesn't affect the kinematics of the universe; that 
is, a universe with Qt = ^o + ^-1/3 has the same H{z) as one with just O^ = Oq (Linder 
1988b).] 

We do not restrict to the flat inflationary case of Oq + ^-1 = 1- We take Tj = 2.73 K, 
i.e. flj = 2.5 X 10~^/i~^, and three neutrino species so O1/3 = 4.2 x 10~^h~^. We ignore 
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Y, take Q^h"^ = 0.0125, and consider h = 0.5 and h = 1 cases. Effects of changing these 
parameters are discussed in Section 5b. 

First, let us explore the early argument that A/ is degenerate with respect to all 
parameters but Or- Initially consider a dust universe. At asymptotically high redshift, 
the angular diameter distance r^ ^ 2H~^ {0)D,~^ z~^ , as can be seen by manipulation of 
(2) into the Mattig (1958; more accessible in Linder 1988b, eqs. 30, A4) analytic form. 
Similarly, from (4) the sound horizon rh -^ (2/^/3)H~^{0)Qo z~^''^ since the baryon 
loading is negligible at very high redshift. Thus A/ -^ tt^/3 z^''^Qo • If one assumes 

— 1/2 

that Zdec = 1100 is in the asymptotic regime then one would obtain A/ ^ I8OO0 
(i.e. ^ ^ 1° X Qo )• For a flat inflationary universe, Oq + fi-i = 1, rh is little affected and 
one asymptotically has Va ^ ^o so A/ ~ O^ - independent of the individual components. 

This sort of estimation led to the early belief that the peak spacing depended only on 
Ot, not Oo and 0_i individually (of course this overlooked the point that by assuming 
a flat universe one already knew Qt ^ what we supposedly were trying to flnd from the 
CMB acoustic anisotropy location). 

However, one is in fact not in the asymptotic regime. In the dust case, the leading 
correction to r^ is of order (0,oz)~^''^, which can be as large as 10%, the influence of 
the radiation component on H{z) is not negligible for rh since Zeq/zdec > (1/20) (it is a 
fairly small effect for Va, < 3%), and the baryon loading effect in Cg is at maximum of 
order {3/4)zeq(b,'y)/zdec ~ 34% (smaller when integrated over z > Zdec)- One actually flnds 
Ta ~ 0~*^'^^ and Th ~ O"'^'^^ so A/ ~ VL~^'^. Remarkably, for a flat inflationary universe 
Oo + ^-1 = 1, J'a ~ ^^°'^°, rh ~ n~°-^^ so A/ ~ 0~°°^ - indeed rather insensitive to the 
individual component densities. This breaks down, however, for nonflat universes. 

In the past few years numerous numerical computations (e.g. Zaldarriaga, Seljak, 
& Bertschinger 1997; Bond, Efstathiou, & Tegmark; Kamionkowski, Spergel, & Sugiyama 
1994), properly including the effects discussed here, have calculated the multipole spectrum 
including the peak spacings. However, no systematic analysis (beyond the useful anima- 
tions of Hu - see fittp://www.sns.ias.edu/ "whu/physics/physics.htmlD has been made of 
the physical scalings, especially for open models, providing simple analytical flts to the 
parameter dependences and demonstrating the role of individual component densities and 
equations of state. Indeed, the myth persists within large sections of the astrophysical 
community that one or two imminent measurements of the flrst few Doppler (acoustic) 
peak locations will flx the total density and hence fate of the universe. Here we attempt to 
portray the results almost as accurately as the matter Boltzmann-radiative transfer nu- 
merical simulations but with explicit flts and scaling of physical dependences, to illustrate 
the true blend of simplicity and complexity in the results. 

Figures 1 show the multipole peak spacing for the various models. Note that the 
vertical spread of the curves demonstrates the lack of degeneracy with respect to the 
individual components - i.e. it is not just Ot that deflnes A/. We do see that flat models 
have a small dispersion, as previously mentioned. The bottom dashed curve gives the 

— 1/2 

naive, often quoted Vtj, ' dependence; this is never an acceptable flt. The bottom solid 
curve gives a derived 0^°'^^ (O^'^'^^ for h = 0.5) flt for pure dust (Og = O^) universes. 
Clearly the presence of individual components does signiflcantly affect the peak spacing. 
Conversely, observations of acoustic peak locations does not uniquely determine the total 
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Figure 1: The peak multipole spacing A/ is plotted against the total cosnio logical density 
Ot for a variety of models. Each panel is labelled by the component equation of state a 
(all include dust as well) and Hubble constant h. The curves from the top down in each 
panel are for constant Oo, starting from Oq = 0.1 (solid), Oo = 0.2 (dotted), Oq = 0.3 
(dashed), etc. They run for values of Ot > ^o- The bottom two curves are fits for pure 
dust models (where Qt = ^o), serving as the lower envelope. The solid bottom curve 
shows the fit of this paper (see equation 7 and Table 1) while the bottom dashed curve 
gives the naive asymptotic fit proportional to Qj, 
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density of the universe. This is discussed further in Section 5a. 

As the equation of state parameter a of the second component approaches zero (that of 
the first component, dust), the dispersion at fixed Ot lessens, as expected. But there is no 



qualitative, smoking gun difference between adding different equation of state components, 
no obvious observational signature. [The pure dust results (Oq = Oj-) are of course the 
same in Figures la and lb, and Ic and Id.] 

Decreasing the Hubble constant has a large effect, increasing both the values and 
dispersion of the peak spacings. Again this is in contradiction to the accepted asymp- 
totic picture, that held that the Hubble constant only entered into the heights and not 
locations of the anisotropy peaks. (Rigorously, A/ would be independent of h only for 
flat universes where Zeq ^ Zdec-) Numerical simulations, however, show the same depen- 
dence exhibited here; for example Figure 2 of Hu & White (1997; also see Figure 3) finds 
A/(0o = 1, ^T = 1)^=0.5 ~ 285, in excellent agreement with our value of 283. Note that 
the dashed asymptotic behavior O^, is here normalized to A/(l, 1) for each h; if we used 
the true naive approach of A/(l, 1) being independent of /i, the fit of this curve to the 
spacings would be even worse. 

Analytic fits to A/ incorporating scaling with the relevant physical variables: the 
matter density Oq, total density O^, and Hubble constant h are useful in presenting the 
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essential physics. Of course in the old flj, fit neither flo nor h would appear. The 
general form we derive is 

A/(o„, Ot)/a/(i, 1) = or + {n-y - or) (jz^) ■ (7) 

Two particular cases of interest are the fiat (possibly inflationary) case when O^ = 1: here 
Al ~ 0~^, and the pure dust case, Oo = Ot, with A/ ~ 0~^. Table 1 gives the values 
X, y, z for the equations of state and Hubble constants adopted. Equation (7) is a central 
result of this paper and provides an excellent flt to the results for A/ plotted in Figures 
1, obtained from solving equations (l)-(5). The approximations are good to about 1% 
rms with maximum deviation about 2%, and thus serve as useful and accurate guides to 
the physical dependences, short of carrying out the full numerical generation of the CMB 
multipole spectrum. 

TABLE 1: Exponents for Peak Spacing Fits 
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-1/3 


0.05 


0.69 


1.18 


can also fit the change 


in peak 


spacing with h as 




A/(0o,0t)/./A/(0„,0t)/^= 


-i = h- 


-/(Ho) 




/(^o) = 


0.30O-°-38, 







(8) 

This is good to better than 2% for h = 0.5. Note that / does not depend on Oy. Also, 
the ratio of the peak spacings for the pure dust case (Oo, Oq) to the fiat case (Oo, 1) are 
independent of /i, depending only on the equation of state, as can be seen from the quantity 
y — X in Table 1. 
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5. Sensitivity to Cosmological Parameters 

We now address how the information resident in the acoustic peaks traces and con- 
strains the cosmological constituents. See Hu & White (1996), especially the excellent and 
comprehensive Section 5.6, for an outline of how the underlying physical process generating 
the matter perturbations influences the anisotropies. 

5a. Determining Cosmological Density 

As illustrated in the previous section, the dependence of the acoustic peak locations 
on the cosmological parameters is more complicated than initially believed. In particular, 
the peak locations and spacing certainly do not completely determine the total density 
Or; there is no unique value derivable without additional assumptions. 

To make this explicit, consider a putative determination of peak spacing at A/ = 350. 
Even if we restrict ourselves to only models with dust and cosmological constant, and 
with h = 1 - i.e. we somehow have prior knowledge that Figure la is the appropri- 
ate one - a horizontal line with A/ = 350 intersects the curves at the following values: 
(no^nr) =(0.48,0.48), (0.4,0.57), (0.3,0.69), (0.2,0.8), (0.1,0.92). This means that a whole 
family of models is allowed, and these have total densities ranging from Qt = 0.48 — 0.92, 
when flo > 0.1. This is very far from a definitive determination of the cosmological density 
and the fate of the universe, even with prior perfect knowledge of the types of cosmological 
constituents and of the Hubble constant. (Even if one restricted to fiat infiationary models, 
note that a measured A/ of 290, say, could arise from Oq = 0.08, 0_i = 0.92, h = 1 or 
Oo = 0.52, 0_i/3 = 0.48, h = 0.5 or O^ = 0.83, 0_i = 0.17, h = 0.5 - quite a range 
of models.) So the next ground based anisotropy detection hinting at a peak should not 
be expected to resolve or even significantly constrain our estimation of the density of the 
universe. 

Lest we be too pessimistic, we note that much additional information exists in the 
amplitudes of the multipole peaks and this will help constrain the cosmological parameters. 
Analysis of these, however, is not so amenable to the quasianalytic arguments employed 
here and best relies on the excellent full numerical codes available. 

5b. Constraining Subsidiary Cosmological Parameters 

Turning from pessimism to optimism, let us consider what we can learn in that bright 
future (when MAP and Planck are aloft?) when superbly detailed CMB and other obser- 
vations have fixed the major cosmological parameters like flo, ^t? h, and a. As mentioned 
in Section 3d a number of other cosmological quantities infiuence the acoustic peak spacing 
in a subsidiary role. The two most promising for constraint are the number of neutrino 
species Ni, and the baryon-photon ratio rj. 

Figure 2 illustrates the sensitivity of A/ to both of these. As discussed in Section 3d, 
the number of neutrino species adds to the relativistic energy density, O1/3, and hence 
affects the redshift of matter-radiation equality, Zeq (which for three species is 21.6 Oq/i^ x 
Zdec)- The other key redshift is that of baryon-photon equality, -2e<j(6,7)5 ^ direct measure 
of rj. The curves are labelled with the baryon loading parameter (3 = {3/4)Zf^q(^f,^^j/zeq = 
0.0158 (O6/iV0.0125) (Oo/i2)-i. 

As one increases the number of neutrino species, A'',^, one decreases Zeq, and increases 
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Figure 2: The peak multipole spacing is plotted relative to what it would be neglecting 
baryon loading and the radiation contribution to the expansion rate. This probes the 
influence of the baryon content of the universe, extra relativistic degrees of freedom, or 
anything that affects the ratio of the redshifts of matter-radiation equality and decoupling. 
The arrows indicate the nominal values (with three light neutrino species) for universes 
with floh'^ = 0.1, 0.25, 1, from left to right. The parameter f3 is proportional to the baryon 
density. Dotted curves neglect the effects on Va to show that the major influence is on 
the sound horizon, dashed curves include the effects on ra for Oq = 1, and solid curves for 



O, 



0.2. 



Al if the baryon loading is small. What is plotted is the ratio of A/ derived with some 
value of Zeq compared to the asymptotic case where the universe at decoupling is completely 
matter dominated {zeq ^ -^dec)- So to find the effect of changing A^,^ one simply takes the 
ratio of the multipole values for the two N^^^s considered. For example, the value at the 
middle arrow (corresponding to three neutrino species in a Qoh'^ = 0.25 dust universe) is 
1.54 for the (3 = 0.01 curve. Adding one light neutrino species decreases log{zeq/zdec) by 
0.055, bringing the value to 1.58, so A/ increases by 2.5%. A reasonable fit is 



{Al)^J{Al)s 



[1 + 0.135 (A^, - 3)] 

0.07(Oo/l2)-0-54 



gi^oh'') 



(9) 



taking into account that (3 changes with Q.oh?. The effects tend to be in the 1-3.5% range. 
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One can also fix the number of neutrino species, and hence O1/3/O0 and z^q/ Zdec^ 
and consider changing the baryon-photon ratio 77, which is proportional to /?. One has 
(3 — [45C(3)/27r^]?B6Xy(Oi/3/Oo) V^ where C, is the Riemann zeta function and rrih = Pb/nb 
the average baryonic particle mass. So if all other parameters were determined, one could 
in principle measure the baryon-photon ratio, or equivalently the specific entropy of the 
universe, s ~ rj~^, from the shift of the acoustic peak spacing relative to its expected value 
for some fiducial 77. Note that even the helium abundance Y enters, through the average 
baryon mass rUb- Ignoring this last, tiny effect, one finds that A/ ~ 7^011 for /3 in the range 
0.01-0.1, leading to a 8% shift for a factor of two difference in ry. 

Recall that for a dust, h = 1 model A/ ~ 0^°'^^ so a 3% change in A/ in this case, 
say, shifts Qt by 5%. The variation of subsidiary cosmological parameters leads only to 
small effects, but the Planck Surveyor should be able to pinpoint the multipole spacing to 
better than half a percent, so they are not wholly beyond detection. 

6. Conclusion 

The physics behind the generation of the CMB acoustic peak anisotropies is well de- 
fined and simple. However, common, asymptotic expressions for the distance behaviors 
and peak spacings are woefully inadequate. Fortunately, there exist comprehensive numer- 
ical solutions, and this paper presents accurate analytic fits to the multipole peak locations 
taking into account the range of relevant cosmological parameters, elucidating the physi- 
cal inputs. While naive application of the peak location to determining the total density 
of the universe is doomed, full use of detailed anisotropy observations does indeed carry 
the hope of deriving not only the total density, but that of the individual components 
- baryons, cosmological constant, etc. - and even such variables as the number of light 
neutrino species and the entropy of the universe. 

This work was supported by NASA grants NAG5-3525, NAG5-3922, and NAG5-4064. 
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